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Abstract 

We develop a formalism to evaluate the Sivers function. The approach is well suited for calcula- 
tions which use constituent quark models to describe the structure of the nucleon. A non-relativistic 
reduction of the scheme is performed and applied to the Isgur-Karl model of hadron structure. The 
results obtained are consistent with a sizable Sivers effect and the signs for the u and d flavor con- 
tributions turn out to be opposite. This pattern is in agreement with the one found analyzing, in 
the same model, the impact parameter dependent generalized parton distributions. The Burkardt 
Sum Rule turns out to be fulfilled to a large extent. We estimate the QCD evolution of our results 
from the momentum scale of the model to the experimental one and obtain reasonable agreement 
with the available data. 
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I. INTRODUCTION 



The partonic structure of transversely polarized nucleons is one of their less known fea- 
tures (for a review, see, e.g., Ref. Nevertheless, experiments for its determination are 
progressing very fast and the relevant experimental effort has motivated a strong theoretical 
activity (for recent developments, see Ref. js}). The present work aims to contribute to this 
effort by using a successful theoretical scenario for the calculation of the Sivers function. 

Semi- inclusive deep inelastic scattering (SIDIS), i.e. the process A(e, e'h)X, with the 
detection in the final state of a produced hadron h in coincidence with the scattered electron 
e', is one of the proposed processes to access the parton distributions (PDs) of transversely 
polarized hadrons. For several years it has been known that SIDIS off a transversely polarized 
target shows azimuthal asymmetries, the so called "single spin asymmetries" (SSAs) 3]. As 
a matter of fact, it is predicted that the number of produced hadrons in a given direction 
or in the opposite one, with respect to the reaction plane, depends on the orientation of the 
transverse spin of a polarized target with respect to the direction of the unpolarized beam. 
It can be shown that the SSA in SIDIS off transverse polarized targets is essentially due 
;o two different physical mechanisms, whose contributions can be technically distinguished 
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7j. One of them is the Collins mechanism, due to parton final state interactions in 
the production of a hadron by a transversely polarized quark [3j| , and will not be discussed 
here. The other is the Sivers is given by 

the product of the unpolarized fragmentation function with the Sivers PD, describing the 
number density of unpolarized quarks in a transversely polarized target. The Sivers function 
is a Transverse Momentum Dependent (TMD) PD; it is a time-reversal odd object [l| and 
for this reason, for several years, it was believed to vanish due to time reversal invariance. 



However, this argument was invalidated by a calculation in a spectator model [9J, following 
the observation of the existence of leading-twist Final State Interactions (FSI) [h]]. The 
current wisdom is that a non-vanishing Sivers function is generated by the gauge link in the 



definition of TMD parton distributions III, Il2|, ll3], whose contribution does not vanish in 



the light-cone gauge, as happens for the standard PD functions. For the same reason it is 
difficult to relate the Sivers Function to the target helicity-flip, impact parameter dependent 
(IPD), generalized parton distribution (GPD) E. Although simple relations between the two 
quantities are found in models Jj], [l^], a clear model independent formal relation is still to 
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be proven, as shown in Ref. 

Recently, the first data of SIDIS off transversely polarized targets have been published, 
for the proton H md toe a— Q. It h as heen foono that, while the S,ve rS effect 
is sizable for the proton, it becomes negligible for the deuteron, so that apparently the 
neutron contribution cancels the proton one, showing a strong flavor dependence of the 
mechanism. Experiments on transversely polarized 3 He target, aimed at extracting the 



neutron information, addressed in 



20. 



2l|. A realistic 



are being performed at JLab 
calculation of nuclear effects for a proper extraction of the neutron information has also been 
performed 22]]. Different parameterizations of the available SIDIS data have been published 
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24 . |25| . still with large uncertainties. Further analyses are in progress (see, i.e. 26|). 



New data, which will reduce the uncertainties on the extracted Sivers function and will help 
discriminate between different theoretical predictions, will be available soon. 

This experimental scenario motivates the formulation of theoretical estimates. One would 
like to perform a calculation from first principles in QCD, however this is not yet possible. 
Lacking this possibility it becomes relevant to perform model calculations of the Sivers 



Q; in the MIT bag 
28 ] and introducing an instanton contribution j^j; in a light- 



function. Several estimates exist, in a quark-diquark model 
model, in its simplest version 
cone model 



301 ]; in a nuclear framework, relevant to establish the manifestation of the Sivers 



function in proton-proton collisions 



3l|. 



To our knowledge, no calculations of the Sivers function have been performed in a Con- 
stituent Quark Model (CQM), i.e. a model described in terms of constituent quarks and 
whose properties have been fixed from hadronic observables. The CQMs have a long history 
of successful predictions in studies of the hadronic spectrum and the low energy electroweak 



structure of 
evolution 



radrons. Ascribing a scale to the model calculations 
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331 ] and using QCD 



351 ] one can evolve the leading twist component of the observable calculated 



in this low energy scale to the high momentum one where DIS experiments are carried 
out. Such procedure has proven successful in describing the gross features of PDs (see, e.g., 
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38]) and GPDs (see, e.g. [39|, |40[), by using different CQMs. Similar expectations 



motivate the present study of the Sivers function. 

In here we propose a formalism to calculate the valence quark contribution to the Sivers 



function from any CQM. Thereafter, we choose the Isgur-Karl model 



41| to perform a 



detailed calculation in order to describe the performance of the approach. A difference 
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in the calculation of TMDs, with respect to calculations of PDs and GPDs, is that the 
leading twist contribution to the one-gluon-exchange (OGE) FSI has to be evaluated. This 
is done through a non-relativistic (NR) reduction of the relevant operator, according to the 



philosophy of constituent quark models 



42|]. 



The paper is structured as follows. In the second section, the main quantities of interest 
are introduced. In the following section, the formalism for the calculation of the Sivers 
function in a CQM is developed. The Isgur-Karl model is presented in the fourth section, 
together with the numerical results of the calculation and their discussion. The following 
section is devoted to the QCD evolution of the model results and to the comparison with 
the available data. In the last section we draw conclusions from our study. 



II. THE THEORETICAL FRAMEWORK 



The Sivers function, f^?{x, kr), the quantity of interest here, is formally defined, ac- 



cording to the Trento convention 
expression 1 : 
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44j, for the quark of flavor Q, through the following 



1 f dj-d\ T 
(2tt) 3 ' 



< T ^(p,s\o Q \p,s) 



(i) 



where St is the transverse spin of the target hadron, the normalization of the covariant spin 
vector is S 2 = —1, M is the target mass and ^(x, kr) is the kr— dependent unpolarized 
PD. The operator Oq is defined as follows 
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13|: 



Q = ^ Q (o, r , 6r)4 (°°> rh + £ (oo, 0)V e (0, 0, 0) 



where iPq(£) is the quark field and the gauge link is: 



(2) 



(3) 



where g is the strong coupling constant. One should notice that this definition for the gauge 
link holds in covariant (non singular) gauges, and in SIDIS processes, since the definition 
of the Sivers function is process dependent. As observed in Ref. [9] for the first time, and 



1 Here and in the following, a ± = (a ± a 3 )/v2 and kr = \kr\ 
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later confirmed using factorization theorems in [45|, |46(, the gauge link, which represents 
the exchange of gluons, provides a scaling contribution which makes the Sivers function non 
vanishing in the Bjorken limit. 

Taking the proton polarized along the y axis one has therefore: 

k T ) = f %^e-^- p+ -^\d Q ) , (4) 



Ak x J (2tt)3 

where the following matrix element has been defined: 

(0 Q ) = { (PSy = l\0 Q \PS y = 1) - (PSy = -l\0 Q \PS y = ~l) } . (5) 

Considering a helicity basis for the target, the Sivers function Eq. (jlj) can be written: 

f&(x,hr) = 9 |^/ d S^ e -^- p+ -^{PS z = l\d Q \PS z = -1) j . (6) 

This equation, finite in the limit of k x — > 0, will be used to evaluate the Sivers function, 
using a CQM to describe the proton. We will now proceed to expand the gauge link, Eq. 
(jSJ), in the coupling constant, g : 

Pexp \-ig A + (r]-^ T )dr]-j = 1 - ig J (77" ,£r)d»T + - (7) 

If the gauge link were not taken into account, it is clear from Eqs. (j2J)- ([6]) that the matrix 
element Eq. (jHJ) would be zero and the Sivers function would vanish. For this reason, the 
first term on the right-hand side of Eq. (CO) does not contribute to the Sivers function. 

A few theoretical predictions have been formulated for the Sivers function. Let us recall 
two of them. 

The first one, based on rather general principles, is the so called Burkardt Sum Rule j^ . 
stating that the total average transverse momentum of the partons in a hadron, (kx), which 
can be defined in terms of the sum of the first moments of the Sivers function for all the 
partons in the target, has to vanish. 

The second one is the conjecture according to which the Sivers function could be related to 
the formalism of the IPD GPDs 48], although, as it has been discussed in the Introduction, 
simple relations between the two quantities are found only in models 



model independent formal relation is still to be proven 



1 1 51 ] and a clear 
161 ] . The IPD GPDs are the Fourier 



transform of the GPDs with respect to the transverse momentum transfer At, at vanishing 
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skewness £. In the case of the helicity independent GPD, Hq(x,£, A 2 ), one has: 



H Q (x,£ = 0,b 2 



d 2 A 7 
(2tt)< 



-ib-An 



H Q (x,£ = 0,A 2 



(8) 



and analogous definitions hold for the helicity independent, target spin-flip GPD Eq(x, £, A 2 ) 
and for the other GPDs. It has been shown that these quantities have a probabilistic 
interpretation, describing the location of the quarks of flavor Q in the transverse plane 
and providing us with a three dimensional picture of the proton {4^. In Ref. 0] 
(see also Ref. 51[ for a recent review on this subject), it has also been shown that, in a 
transversely polarized proton, for example along the y direction, the quantity describing the 
distribution of the quarks of flavor Q, with longitudinal momentum x, in the transverse 
plane, independently of their helicity, is 



1 b S d 

Pq(x, e = 0, b) = -H Q (x, 0, b 2 ) - -^JL—E Q (x, 0, b 2 ) 



(9) 



i.e., the transverse polarization of the proton produces a shift in the transverse location of 
the quarks. As explained before, this effect in the partonic structure of transversely polarized 
protons has been related, in peculiar models, in a qualitative way, to a nonvanishing Sivers 



effect 



49 



50]. 



III. THE SIVERS FUNCTION IN CONSTITUENT QUARK MODELS 



The constituent quark, one of the most fruitful concepts in 20th century physics, was 
proposed to explain the structure of the large number of baryons being discovered in the 



sixties 52J. The constituent quark concept was incorporated into a QCD scheme by taking 
into account gluon exchanges 421 ] . The chosen description was a potential model in order 
to establish an immediate connection with all previous work. 

The constituent quark scheme has guided some of the most successful parameterizations 
of parton distributions |53|. Besides, the philosophy that has guided these parameteriza- 
tions is precisely the one used to establish the link between constituent models and parton 
distributions. More specifica 



ly, model calculations are ascribed to a scale determined by 
their partonic content 32j, |33(. In most models that scale is characterized by the existence 
of valence quarks only. From that low scale one uses DGLAP evolution to describe the 
partonic regime 



The models based on constituent quarks (CQMs) have produced beautiful results in the 



description of PDs and GPDs, leading to a phenomeno 
terms of momentum densities and wave functions 
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ogical understanding of them in 
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401 ] . This success in the 



description of many parton distributions makes us confident that the application of the 
approach to the calculation of the Sivers function will also serve to guide the experimental 
observations and help the physical interpretation of this observable. 

Let us specify in detail the scheme in which we are going to develop our formalism for the 
Sivers function. We shall assume that the nucleon at a certain low energy scale is made up 
of valence quarks only. These valence quarks are held together by a confining interaction; in 
addition, there is a residual interaction, governed by the structure of perturbative QCD, e.g. 
the One Gluon Exchange Interaction. The strong confining interaction maintains the quarks 
together, while the residual one governs the splittings within the same flavor multiplet. Any 
scheme with these hypothesis is a constituent quark model framework. 

This scheme has never to be understood in a trivial perturbative sense. The parameters 
absorb much of the non perturbative features of the dynamics and this relation between the 
parameters and some chosen observables makes the scheme predictive. If one goes to higher 
order in the perturbative expansion, one needs to find new parameters from the chosen 
observables. Thereafter, the predictions do not change much with respect to the lowest 
order result [54]. Certainly we are dealing with models and not with QCD and therefore one 
should not expect precision. Nevertheless, the scheme has been so successful that particles 
which do not fit approximately under it are called exotics, hybrids or other peculiar names. 

Using this scheme we evaluate a formula for the Sivers function, defined according to 
Eq. ([6]), valid for any CQM. Let us proceed to the analysis having in mind Fig. 1. To the 
first non vanishing order giving a contribution to the asymmetry, the Sivers function for the 
flavor Q is obtained as follows: 

fjf(x, k T ) = 9 { iL f ^f^e-^-P^^kr)^ ) ; (1Q) 

where 



(d Q ) = (p^ = i|^(o,e-,&)(^){o a (o,r,&)^} 

x 7 + ^(0)|PS z = -l)+h.c. , (11) 
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where T£j = A^-/2 with A?- being a Gell-Mann matrix, and 

_, POO _, 

0«(O,£- 6r) = / A+(0, V -,fr)d V - 



In the above equations, use is made of light-cone states 2 , defined as \p) = \p + ,pr), with 
p- = (m 2 +p 2 T )/(2p + ). The light -cone states are normalized as follows: 

(p'r'lpr) = (2n) 3 2p + 5(p' + - p + )5{$ T - p T )8 rr > , (13) 

where the label r represents a set of discrete quantum numbers. The creation and annihila- 
tion operators of the quark fields are normalized accordingly: 

{6f(p), h(p')} = (27T) 3 2p + 5(p' + - p + )5{p T - p T )5 lv , (14) 

where the set I = {m, c, J 7 } includes the helicity, color and flavor quantum numbers of the 
quark, respectively. 



Using the approximation of expanding Eq. (iTTj) in terms of free quark fields |39|] , one gets 

X / 4Efi(fc3)^ r - lfc ' rT «m 3 (4) 
J m 3 

x ^){6 a (o,r,fT)^}7 + 

x £/ d&b^kQu^ii&lPrS, = -1) + h.c. 1 , (15) 

where dk t = dkfdkTi/{2kf(2i[) 3 ). Inserting now proper complete sets of intermediate free 
one quark states, the previous equation becomes 



b 3j 

ma 

x (w) £ / d ^n j dk^kj^lkjn) (k n l n \(kili\ 



2 Here and in the following, x = (x + ,x ,xt) is a four vector in light-cone coordinates, while obviously 
x = (xi,X2,Xa) and x T = (xi, x 2 )- 
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x {6 a (o,r,&)^} 

x E / dk' n f d^|^ii>|^X>fel(^^l7 + 

x E/ = -l) + h.c.\ . (16) 
If there is no further interaction within the recoiling system, one has: 

(kJnKO = (2rr) 3 2k+6(k> n + - k+)5(k n , T - k nT )5 lnA , (17) 

(PrS z = l | {b%{h)\hh)\k n ln)} 

= (2ixf2k + J(P + - kf - k+ - K)5{P T - k 1T - k 3T - k nT )5 {Sz , rMMM 
x (PrS z = 1| k 3 {m 3 ,i, Q}; k 1 {m 1 ,c 1 ,J r 1 }; P - k 3 - hjn) 
= (27t) 3 2k+5(P + - k\ - k+ - k+)5(P T - k 1T - k 3T - k nT )5 {Sz ,r,i M 
x ^t Sz =i (k 3 {m 3 , i, Q}; k 1 {m 1 , c Xy P — k 3 - k%, Z n ) . (18) 

In the last equation, the definition of the intrinsic proton wave function, in momentum 
space 3 , has been recovered. In the same equation, the terms Srs z ,r,...) are showing that all 
the discrete quantum numbers of the quarks have to be properly combined to recover those 
of the parent proton. In order to obtain a workable expression for the Sivers function given 
by Eq. (TT6l) . other three relations have to be used. One is written using Eq. (|T2|) and 
translational invariance: 

(Mil {d a (0,C,£ T )} WA) = e^^-^^Ck.kl {6 a (0)} \k[l[)e- ik '^- +i ^ . (19) 

is the identity Q = 

0.(0) = f AU r ,0 T ) ir = - / ^^-AM . (20) 

The last one is obtained by evaluating the matrix element of the perturbative free gluon 
operator appearing in Eq. ffl~9]) . Assuming, as an approximation, that this operator is 
time-independent, one gets, in the Landau gauge: 

x (2Tif2k+5(k+ - k'+ - q + )5(k 1T - k' 1T - q) . (21) 

3 In the class of models to be later used, the separation of the center of mass and intrinsic motion is always 
possible. 
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Substituting in Eq. ([TBI) the identity: 

1 1 



. = i(2n)5(q + ) , (22) 

together with Eqs. (JTTj) - (I2TI) . one is left with the following expression for the Sivers 
function: 

f^f(x,k T ) = ^ j-^ 2 7^- / d kidh-j7^S(q + ) 

x + q + _ a;P + )5(fc3T + <?t - M(27r)5(go) 

JFi ,mi ,c\,m!^ ,Cj ,m3,i,m' 3 ,j 

x * r Sz =-i (fcs + g, {^3. i> Q}; fa - & Wi> c i) ^i}; P-fo- h, k) } ,(23) 

with the interaction determined by: 



V(fci, fc 3 , g) = -^M m3 (4)7 + w m ' (h + g) u mi(^i)7 + «m' (&i - <f) ■ 
qz 3 i 



(24) 



Since the final aim is the evaluation of the Sivers function within a NR model, a NR 
reduction of the interaction has to be performed. Using therefore the definitions of free 
four-spinors in Eq. ( 1241) . performing a NR expansion leaving out terms of second order in 
momentum, as it is commonly done in nuclear physics (cf. Ref. 55|]), one gets the potential: 

1 



V NR (ki,k 3 ,q) 



2q 2 



(Vo) mi,m. 



m3,m 3 



(25) 



with: 

Vo(ki, k 3 , q)m 1 ,m' 1 ,m 3 ,r. 



kl t q-k 3 i k\ q- k x 
4m 2 



1 + ^ + ^ + ^--^ + ° 

m Am z 



ni 









m 2 J 



&mi,rn l 1 &m3,m', i (26) 



Vs{ki, k 3 , gOmi.mi, 



-;|1 + S + £^ M . 
m Am z 



[q X (gQmi.m'Ji 

2m 



kf_ q- kA [q x (a 3 ) m3 ,m' 3 ]z 

m 4m^ / 1 2m 
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1 4 m 2 ^m 3 ,m^ 

[g x (ai) ■ (fc 3 x q) 
_l i ^ 

8m 3 

mi ,m, 

■ (fci x x (<t 3 ) 

ma ,m Q I z 

_| i i — 

8m 3 



4m 2 



A few remarks are in order. First of all, the helicity conserving part, Vq, Eq. fl26|) . of the 
global interaction Eq. (I25I) . does not contribute to the Sivers function. One should notice 
that, in an extreme NR limit, the Sivers function would turn out to be identically zero. In 
our approach, it is precisely the interference of the small and large components in the four- 
spinors of the free quark states which leads to a non-vanishing Sivers function, even from 
the component with I = of the target wave function. Effectively, these interference terms 
in the interaction are the ones that, in other approaches, arise due to the wave function (see, 



e.g., the MIT bag model calculation 



28|) 



The scheme is now completely set up and any CQM can be used to evaluate the Sivers 
function. We next use properly normalized NR wave functions to transform Eq. ( |23i) in: 

f£?(x, hr) = 3? "^ 2 — J dhdh-^-Sikt - xP + )5(k 3T + q T - k T )M Q , (28) 



where 



M Q = £ ^^x^^,^,;,^) 

T\ ,mi ,ci ,m'j ,Cj ,m3 ,i,m' 3 ,j 

x ^ts z =i (k 3 {m 3 ,i, Q}; ki{m u c u P-k 3 - k h l r . 
x T^V^k^q) 

x ^ rSz =-i (k 3 + q, {m' 3 ,j, Q}; h - q, {mi, c'^Fi}; P - k 3 - kx, l n ) . (29) 

Each wave function ^ r s z describing a possible proton state can be factorized into a com- 
pletely antisymmetric color wave function, x, and a symmetric spin-flavor-momentum state, 
as follows: 

= ®sf,s z (k 3 {m 3 , Q}; fcihi, T x };P-k 3 - k u {m nj x(h Ci, c n ) . (30) 

11 



The matrix element of the color operator in Eq. (1291) can be therefore immediately evaluated: 

£ X \h ci, c^T^xU, c[, c n ) = -| , (31) 

ci,c[,i,j 

which is the well-known result for the exchange of one gluon between quarks in a color 



singlet 3-quark state [56( . Besides, as a consequence of the symmetry of the state one 
can assume that the interacting quark is the one labeled "3", so that, after the evaluation 
of the summation on the flavors T\, M. a can be written, for the u and d flavors, as follows: 

M u(d) = . 3 . ^ $l f ,s.=i (£ 3) mjj . £^ mi . p _ ^ _ ^ mn 

m\ ,771^,7713,7773 



x l±r 3 3 y NR (J^ u %^^ 

x $ s/jSz= _ 1 (& 3 + q, m' 3 ; k\ - q,m[; P -k 3 - h, m n ) . (32) 

Eq. (128]) . with A4 U ^ given by Eq. (1321 . provides us with a suitable formula to evaluate the 
Sivers function, once the spin-flavor wave function of the proton in momentum space, i.e. 
the quantity $ s j, is available in a given constituent quark model. 

IV. THE CALCULATION OF THE SIVERS FUNCTION IN THE ISGUR-KARL 
MODEL 

As an illustration, in this section we present the results of our approach in the CQM of 



Isgur and Karl (IK) |4lj. In this model the proton wave function is obtained in a OGE 
potential added to a confining harmonic oscillator (H.O.); including contributions up to the 
2hu shell, the proton state is given by the following admixture of states 

\N) = a\ 2 S 1/2 ) s + b\ 2 S[ /2 ) s + c\ 2 S l/2 ) M + d\ A D l/2 ) M , (33) 

where the spectroscopic notation \ 2S+1 Xj) t , with t = A, M, S being the symmetry type, has 
been used. The coefficients were determined by spectroscopic properties to be a = 0.931, 



b = -0.274, c = -0.233, d = -0.067 [56j. If a = 1 and b = c = d = 0, a simple H.O. model 
is recovered. The parameter a 2 = moo of the H.O potential is fixed to the value 1.23 fm~ 2 , 
in order to reproduce the slope of the proton charge form factor at zero momentum transfer 



56|. 
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The formal expressions of the wave functions appearing in Eq. (1331) in the IK model can 



be found in 



56, 



571 ]. given in terms of the following sets of conjugated intrinsic coordinates 



R = ^(fi + r 2 + r 3 ) 



P 



V2 



(rl - r 2 ) 



A 



(rl + f 2 - 2r 3 ) 



K 

k p 
kx 



-^(h + h + h) , 



V2 
1 



(ki - k 2 ) , 
(kt + k 2 - 2k 3 



(34) 



There are many good reasons to use the IK model as a test of the developed formalism. 
First of all, the IK is the typical CQM, succesful in reproducing the low-energy properties 



of the nucleon, such as t 
momentum transfer 41 



re spectrum and the elastic and transition form factors at small 



561 ] . In particular, as was shown in Ref. [581 ]. in the IK model, 



(k 2 )/m 2 ~ 0.3 and therefore one expects small corrections from terms O (k 2 /m 2 ). Besides, 



one of the features of the IK model is that the OGE mechanism [42j, which reduces the 
degeneracy of the spectrum, is taken into account. It is therefore natural to study our 
formalism, based on OGE FSI, within the IK framework. Concerning PDs, it has been 
shown that the IK model can describe their gross features, once QCD evolution of the 
proper matrix elements of the corresponding twist-2 operators is performed from the scale 



of the model to the experimental one 



36 



37 



. Reasonable predictions of GPDs have 



also been obtained [39], and this makes particularly interesting the evaluation of the Sivers 
function in the IK model. In section II, the relation between the Sivers function and the 
impact parameter dependent GPDs has been discussed. In a model where a shift of the 
quark location in the transverse plane is found, a sizable Sivers function should arise. In 
order to investigate whether the IK model is suitable for the analysis of the Sivers function, 



performing the Fourier 
39j. The quantity: 



the quantity Pq(x, £ = 0,6) has been calculated in this model 
transforms, Eq. ([8]), of GPDs evaluated along the lines of Ref. 

PqQ>) = J dxp Q (x,£ = 0,6) , (35) 

representing the distribution of the quarks of flavor Q, with any longitudinal momentum, in 
the transverse plane, independently of their helicity, in a proton polarized along the positive 
y direction, is shown in Fig. 2. It is clear that a slight shift along the x direction is observed, 
with a different sign for the u and d flavor. Therefore, according to the present wisdom, a 
small Sivers function is expected, with different sign for the u and d flavors [49J . 
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After this discussion, the IK model appears as a promising framework for the evaluation 
of the Sivers function. 

The Sivers function has been calculated according to Eq. (1251) . using the proton states 
Eq. (1551) . neglecting the small D component, and the potential Eq. (1231) in A4 U ^ given by 
Eq. dS2D. 

The results of the calculation can be cast in the following form: 

Q=u4 . 1 . V2g 2 M 2 (3\l 1 f d 2 q T k° x -£fe+S4--s/!**A 



f 1T ' (x, hr) - ^ [ 2 ) 27r3/2a 3 J (2?r)2 |fc o _ e 



x 



1 

2m 



(36) 



with k® = Jm? + k\, and 



x 



V2 W ; ' 2V3P+X 



_ k z \ + ^at • (37) 

The expressions of the functions p*xx are given in the Appendix. 

To evaluate numerically Eq. (1361) . the strong coupling constant g, and therefore a s (Q 2 ), 
has to be fixed. Here, the prescription introduced in the past for calculations of PDFs in 



quark models (see, i.e., Ref. 38j) will be used. It consists in fixing the momentum scale of 
the model, the so-called hadronic scale according to the amount of momentum carried 
by the valence quarks in the model. In the approach under scrutiny, only valence quarks 
contribute. Assuming that all the gluons and sea pairs in the proton are produced pertur- 
batively according to NLO evolution equations, in order to have ~ 55% of the momentum 
carried by the valence quarks at a scale of 0.34 GeV 2 , as in typical low-energy parameteriza- 
iions [53J, one finds, that \i\ ~ 0.1 GeV 2 if Ag£g ~ 0.24 GeV. This yields a s (/i 2 )/(4vr) ~ 0.13 



38]. 



For an easy presentation, the quantity which is usually shown for the results of calculations 
or for data of the Sivers function is its first moment, defined as follows : 

-(1)2/ n _ f Jit k T 



ft^(x) = I d 2 k T ^-J^(x, hr) . (38) 
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The results of the present approach for the moments Eq. (1381) are given by the dashed 
curves in Fig. 3 (4) for the u (d) flavor. They are compared with a parameterization of 



the HERMES data, corresponding to an experimental scale of Q 2 = 2.5 GeV 2 24j 4 The 



24] 



patterned area represents the 1 — a range of the best fit proposed in Ref. 

As expected from the IPD GPDs analysis, shown in Fig. [2J a different sign for the u and 
d flavor is found. 

Let us see now how the results of the calculation compare with the Burkardt sum rule 



471 ]. which follows from general principles and must be satisfied at any scale. If the proton 
is polarized in the positive y direction, in our case, where only valence quarks are present, 
the Burkardt sum rule reads: 

E <^ Q > = , (39) 

Q=u,d 

where 

<*?> = - £ dx J dk T ^tf T Q (x, k T ) . (40) 

Within our scheme, at the scale of the model, it is found = 10.85 MeV, (k^) = 

— 11.25 MeV and, in order to have an estimate of the quality of the agreement of our results 
with the sum rule, we define the ratio 

_ (K) + (K) / 41 \ 
(k*) - (kf) ' 1 ] 

obtaining r ~ 0.02, so that we can say that our calculation fulfills the Burkardt sum rule to 
a precision of a few percent. 

Another prediction has been derived in the framework of large N c (goJ and it reads, when 
xN c ~ 0(1) and the large N c predictions are supposed to be applicable: 

We get the closest value to the prediction above, 0.26, in a narrow region around x = 0.4. 



It has been chosen to compare the results with the parameterization of 24| and not with that of 23| or 



251 ] just because, in the first case, it is easier to reconstruct the parameterization of the data, and their 
1-sigma range has been kindly provided by the authors of Ref. [24] . The discussion of the quality of the 
agreement of the present results with data would not change substantially if the comparison were made 



with the parameterization of Refs. 23|, [25] 
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We note that the contribution of the states \ 2 S'y 2 )s an d | 2 «Si/2)m; m spite of their small 
probability in the proton state Eq. fl33l) . turns out to be important in the evaluation of the 
Sivers function. 

The magnitude of the results is close to that of the data, although they have a different 
shape: the maximum (minimum) is predicted at larger values of x. One should anyway 
realize that one step of the analysis is still missing: the scale of the model, /j,q, is much lower 
than the one of the data, which is Q 2 = 2.5 GeV 2 . For a proper comparison, the QCD 
evolution from the model scale to the experimental one would be necessary. This issue is 
discussed in the next session. 



V. QCD EVOLUTION OF THE MODEL CALCULATION 

The Sivers function is a TMD PDs and the evolution of this class of functions is, to a 
large extent, still to be understood. In any case, recent interesting developements can be 
found in Ref. (m). 

In order to have an indication of the effect of the evolution, we perform a NLO evolution 
of the model results assuming, for the moments of the Sivers function, the ones defined 
in Eq. (138]) . the same anomalous dimensions of the unpolarized PDFs. As described in 
the previous section, the parameters of the evolution have been fixed in order to have a 
fraction ~ 0.55 of the momentum carried by the valence quarks at 0.34 GeV 2 , as in typical 



parameterizations of PDFs j53|, starting from a scale of /Aq — 0.1 GeV 2 with only valence 
quarks. The final result is given by the full curve in Fig. 3 (4) for the u (d) flavor. As it is 
clearly seen, the agreement with data improves dramatically and their trend is reasonably 
reproduced at least for x > 0.2. 

Of course a word of caution is in order: the performed evolution is not really correct. In 
any case, an indication of two very important things is obtained: 

i) The evolution of the model result is necessary to estimate the quantities at the mo- 
mentum scale of experiments, as it happens for standard PDs 36l. 1371. l38l|; 

ii) after evolution, the present calculation could be consistent with data, at least with the 
present ones, still affected by large statistical and systematic errors. 
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VI. CONCLUSIONS 



A rather general formalism for the evaluation of the Sivers function, to be used in any 
CQM, has been developed. The crucial ingredient has been the NR reduction of the leading 
twist part of the OGE diagram in the final state. It has been shown that the IK model, based 
also on a OGE contribution to the Hamiltonian, is a proper framework for the estimate of 
the Sivers function. The obtained results show a sizable effect, with an opposite sign for 
the u and d flavors. This is in agreement with the pattern found from an analysis of impact 
parameter dependent GPDs in the IK model. 

Let us compare our approach with previous calculations. The diquark model with scalar 
diquarks has no contribution for the <i-quark [^tJ and therefore does not satisfy the Burkardt 
sum rule (BSR), Eq. fl39l) . The diquark model with axial- vector diquarks has contributions 
to both u and <i-quarks and with opposite sign, but with the magnitude of the d 10 times 



smaller than that of the u. The BSR is not satisfied. The MIT bag model calculation [28] 

has non-vanishing u and cZ-quarks contribution of opposite sign which are proportional in 

magnitude. The <i-quark contribution is much smaller than ours and therefore does not 

[ ] 

satisfy the BSR. The MIT bag model modified by instanton effects ]29] has u and <i-quark 
contributions of the same sign and therefore does not satisfy the BSR. As a summary, we 
can say that our calculation, despite the naive wave function used, is in better agreement 
with the data with respect to the other approaches, and fulfills the BSR. 

In order to compare with the data, one has to evolve the model calculation to the experi- 
mental scale. Although a consistent QCD evolution of the model results to the experimental 
momentum scale is not yet possible, due to the lack of the calculation of the corresponding 
anomalous dimensions, an estimate of the evolution has been attempted. It has been found 
that, once properly evolved, the model results could be in reasonable agreement with the 
available data. 

The formalism presented here can be used with any CQM and it will be interesting in 
the near future to implement other calculations with different models, performing a correct 
evolution as soon as the corresponding ingredients become available. The connection of the 
Sivers function with IPD GPDs deserves a careful analysis and will be discussed elsewhere. 
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APPENDIX: THE SIVERS FUNCTION IN THE IK MODEL 



The functions appearing in Eq. (|36|) are listed below. 
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FIGURE CAPTIONS 



Fig. 1: The contributions to the Sivers function in the present approach. The graph has 
been drawn using JaxoDraw 62]. 



Fig. 2: In the upper (lower) panel, the quantity Pq(&), Eq. (I35I) . is shown for the u (d) 
flavor. 

Fig. 3: The quantity fyr ( x )i Eq. ( |38l) . for the u flavor. The dashed curve is the result 
of the present approach at the hadronic scale /xjj, Eq. (1361) . The full curve represents the 
evolved distribution after standard NLO evolution (see text). The patterned area represents 



the l — o" range of the best fit of the HERMES data proposed in Ref. 



24]. 



Fig. 4: The same as in Fig. 3, but for the d flavor. 
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